Universal mechanism of dissipation in Fermi superfluids at ultra low temperatures 
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We show that the vortex dynamics in Fermi superfluids at ultra-low temperatures is governed by 
the local heating of the vortex cores creating the heat flux carried by non-equilibrium quasiparticles 
emitted by moving vortices. This mechanism provides a universal zero temperature limit of dissipa- 
tion in Fermi superfluids. For the typical experimental conditions realized by the turbulent motion 
of 3 He-B the temperature of vortex cores is estimated to be of the order 0.2T C . The dispersion of 
Kelvin waves is derived and the heat flow generated by Kelvin cascade is shown to have the value 
close to the experimentally observed. 

PACS numbers: 



Vortex dynamics determines the most fundamental 
properties of supcfluids and superconductors. The well 
known textbook example is the motion of Abrikosov vor- 
tices under the action of the Lorentz force which deter- 
mine the finite resistance of the type-II superconductors 1 -. 
Equally important is the vortex dynamics in the super- 
fluids 4 He and 3 He where motion of quantized vortices 
mediate the mutual friction force between the normal 
and superfluid components originating from the scatter- 
ing of normal excitations by moving vortex lines. The 
dissipative component of mutual friction provides the re- 
laxation of superflow which is the key aspect in the the- 
ory of superfluid turbulence^. Recently there has been a 
renewal interest in this field owning to the developing ex- 
perimental techniques allowing to study in particular the 
decay of quantum turbulence at very low temperatures 3 -. 
In this regime the normal component of the fluid be- 
comes extinct and the mutual friction can not provide 
the energy dissipation and the relaxation of the super- 
flow. The state-of-art experiments 4 - - — demonstrate that 
the dissipation in the vortex motion remains finite at the 
lowest temperatures both in the Bose superfluid 4 He and 
Fermi one 3 He-B. The energy dissipated by turbulence 
is supposed to be released in the form of the thermal 
flux of non-equilibrium quasiparticles which recently has 
become accessible to direct measurements^. These ex- 
perimental results pose a challenging question about the 
fundamental nature of the dissipation in superfluids at 
ultra low and even zero temperatures. 

In this Letter we propose the new mechanism which 
governs the vortex dynamics in Fermi superfluids at 
ultra-low temperatures. We show that even at zero tem- 
perature in the absence of the normal component the 
non-stationary dynamics of vortices is intrinsically dissi- 
pative as a result of relaxation processes inside the vor- 
tex core moving with the finite acceleration (which is a 
generic situation in superfluid turbulence). We address 
the case of typical superconductors and superfluid Fermi 
systems which can be described within the weak coupling 
Bardeen-Cooper-Shcriffer (BCS) theory. In particular 
such restrictions means that the critical temperature and 
energy gap are much less than the Fermi energy which 



allows the so-called quasiclassical description of the pair- 
ing stated. Note that for the case of cold atomic gases 
this description is valid only in deep BCS limit. 
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FIG. 1: (a) Container with the ideal gas of particles demon- 
strating the classical analog of vortex core. The temperature 
of gas is raised in case if u / 0. The energy dissipation 
by the heat flux through the heat conducting wall is shown 
schematically, (b) Sketch of the vortex core and quasipar- 
ticle confinement due to Andreev reflection and the flow of 
non-equilibrium delocalized quasiparticles. 

The mechanism of the ultra-low temperature energy 
dissipation in Fermi superfluids has a very transparent 
classical physics analogy. As will be discussed below the 
vortex core can be viewed as an insulating box containing 
a gas of monoatomic molecules which represent localized 
vortex core quasiparticles shown in FigQJa). Let us as- 
sume that the oscillating center of mass motion of the box 
is produced characterized by the velocity u(t) = u Q(t), 
where Q{t) is a step function with a period t p larger 
than the thermalization time of the gas. Then each step 
change of the velocity will result in the temperature rise 
which can be calculated from the energy conservation law 
to be ATi oc = 4uq/c v where c v is the specific heat per 
unit mass. Averaging over the time interval much larger 
than the period we obtain the linear temperature growth 

in time TJ oc (t) = -^-j-- Furthermore, the temperature 
growth will certainly stop at some point regulated by 
the finite heat conductivity of the container walls. The 



2 



balance of the heating by box acceleration and the heat 
release through the heat conducting walls will set the 
stationary value of the temperature T/ oc of the particles 
localized inside the box. The heat flow through the con- 
tainer walls determines the rate of the energy dissipation 
in the system. 

At temperatures considerably lower than T c vortex dy- 
namics in clean type-II superconductors and Fermi su- 
perfluids is determined by the kinetics of localized ex- 
citations bound in the vortex cores driven out of equi- 
librium by vortex motion^. In Fermi superfluid contain- 
ing vortices the profile of the order parameter A (r) near 
the vortex core produces a potential well where localized 
states with a discrete spectrum exist. The quasiparticle 
confinement results from the subsequent Andreev reflec- 
tions transforming particles to the holes and vice versa 
[FigQJb)]. The localized states correspond to energies 
|e; oc | < Ao where Ao is bulk value of the gap. The spec- 
trum has the so called anomalous branch^ whose energy 
varies from — Ao to +Ao and crosses e = 0. The anoma- 
lous energy branch has topological origin^ and exists 
for all kinds of quantized vortices in the weak coupling 
Fermi superfluids which allows for the quasiclassical de- 
scription. In model axially symmetric singly-quantized 
vortices the anomalous branch states are characterized 
by the particle angular momentum projection on the vor- 
tex axis /i. For low energies |e/ oc | <C Ao, the anomalous 
branch is e; oc = —tkuofi. The separation between the 
levels with neighboring angular momenta /j, and /U + 1 is 
luq ~ Ao /(hkp£) where kp is the Fermi momentum and 
£ is the coherence length. In the weak coupling limit the 
interlevel distance is small ljq <C Ao therefore the angular 
momentum [i can be considered as classical variable com- 
muting with the corresponding angle 9 P , characterizing 
the direction of quasiparticle motion in plane perpendic- 
ular to the vortex axis. 

Our results remains valid with some modifications for 
the non-axisymmctric vortices existing in multicimpo- 
ncnt Fermi superfluids like 3 HoAS.. In this case one 
should assume the angular dependence of bound states 
energy a oc = e; oc (^, fl p ) n ' 13 . Note also that in general 
the singly-quantized vortices have two anomalous energy 
branches corresponding to different spin states of bound 
fcrmions— . The kinetic theory that we will use allows to 
treat each of the anomalous branches in separate. There- 
fore below we omit the spin indices, which finally reduces 
the full Bogoliubov-dc Gcnncs 4x4 matrix theory to the 
effective scalar theory of the core fermions. 

The key idea of the present work is that the ensemble of 
quasiparticlcs localized within the moving vortex core is 
in some sense analogous to the particle gas confined in a 
moving box. Below we will show that the local heating of 
the localized quasiparticlcs due to the vortex acceleration 
can lead to the temperature growth which is compensated 
by the heat flow out of the vortex cores in the form of 
the non-equilibrium delocalized quasiparticlcs. 

The equilibrium distribution function of Fermi quasi- 
particles has the form /(°) = tanh(e/2T). The delocal- 



ized quasiparticlcs with energies |e| > Ao have the tem- 
perature T = Td e i fixed by the heat bath. We assume 
that Tdei — > which allows to neglect the concentration 
of delocalized particles so that /<S(e > Ao) = 1 and 

fdei( e < ~ A-o) = — 1. On the other hand as we will see 
below the localized particles have the different tempera- 
ture T = Ti oc which can be much larger T\ oc 3> Tdei ■ 

We start with a kinetic equation^ for the distribution 
function of localized quasiparticlcs / = f(t,/j,,9 p ) which 
depends on time t and a pair of canonical variables (fj,, p ) 

_„ o([pxu] . z) | + | + v ej+ g A=S((/) (1) 

Here the energy spectrum of localized particles e; oc = 
— ftwoM plays the role of the classical Hamiltonian gov- 
erning the dynamics of the canonical variables so that 
8 p = —luq and fi = 0. The source of non-equilibrium is 
the first term in Eq. ([I} which is generated by the vortex 
motion with the velocity u = (u x x, U v y). 

The collision integral (CI) in the r.h.s. of this equation 
(pTJ) can be determined by different scattering processes 
which depend on the particular system where the Fermi 
superfluidity is considered. We will consider the universal 
mechanism which inherently exists in any Fermi liquid 
- the mutual scattering of quasiparticlcs. In 3 He and 
ultra-cold Fermi gases the collision integral is determined 
ultimately by this mechanism. 

The general form of the particle-particle CI has rather 
complicated form-^. Therefore we use the model forms 
of the CI which features however the important proper- 
ties of the exact CI. We are interested in the relaxation 
of localized quasiparticle with the energy |e| < Ao and 
momentum p. It collides with a particle characterized by 
ei, pi and they scatter into states 62, P2 and £3, P3. Dur- 
ing this process the energy and momentum are conserved 
e + ei = e 2 + e 3 and p+pi = P2 + P3- Due to the Gallilean 
invariance the CI integral vanishes when all quasiparticlcs 
are moving with a net velocity w having the equilibrium 
distribution function = tanh[(e — pw)/2T]. 

In general, there are three types of scattering processes 
which contribute to the particle- particle CI which we 
consider in separate. 

(i) Collision of localized and delocalized qusiparticlcs 
when |e|, |ei| < Ao and at the same time | £2 .3 1 > Ao. 
This is the basic scattering process determining the force 
acting on the vortex from the heat bathi&. For T^ei — > 
the probability of this collision is negligibly small and for 
Tdei = the exact calculation shows that the correspond- 
ing part of the CI is zero. 

(ii) The scattering involving only localized quasiparti- 
cles |e|, I e 1 1 , | £2 1 , £3 < Ao. In this case the collision time 
approximation can be used with account of the Gallilean 
invariance in the system of localized particles^ 



St loc (f) = — [/(°)( e -pw)-/]. 

Tloc 



(2) 



Here w is the 'drag velocity' which determines the evolu- 
tion of distribution function. In the absence of collisions 
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with dclocalizcd states it coincides with the vortex veloc- 
ity w = u. The collision time depends on the temper- 
ature of localized quasiparticles and for T\ oc <C A was 
estimated to be±& r^ 1 — t~ 1 (Ti oc / A ) 4 where r^ 1 — 
Aq/Ef is the normal state scattering rate at T = T c . 
Importantly this part of CI involving only localized par- 
ticles © conserves the energy of localized particles 



A 



eSt loc (f)de = 0. 



(iii) Collision of two localized quasiparticles with pro- 
duction of localized and delocalized ones. In this case 
we assume that |e|, |ei|, |e 2 | < A and |e3 1 > A . During 
such processes the non-equilibrium delocalized quasipar- 
ticles are created which carry the energy flow out of the 
vortex cores shown schematically in Fig 0Tb). This mech- 
anism of energy losses is similar to the heat flow through 
the heat conducting wall in the gas container discussed 
above [see FigQJa)]. 

To calculate the probability of scattering event (iii) 
one needs to take into account the contributions from 
states with £3 > Ao and £3 < — A<j. We will con- 
sider in detail only the first case and the second one 
can be treated analogously. We can put f^°'(ea) = 1 
which means that the final delocalized state is empty. 
The scattering probability is determined by the factor 
F(e,e 1: e 2 ) = [1 -/(e)][l - /(ei)][l +/(e 2 )j which is the 
product of the initial states occupational numbers and 
the probability of the final localized state to be empty. 
To obtain the total scattering probability we should in- 
tegrate over the energies taking into account the en- 
ergy conservation requirement — Ao < £2 < e* where 
e* = E\ +e — Ao which yields in particular Ao > £1 > — e. 
Then we obtain the following contribution to CI 



Sti- d i(f) = (r, 



nA^y 1 [ ° da [ de 2 F(e,e 1 , 

J-e J-Ao 



e 2 ) (4) 



where the prcfactor can be determined from exact form 
of the CI. The contribution from the states with £3 < 
— Ao yields the same result Eq.© up to the numerical 
prefactor which depends on the particle-hall asymmetry. 
Note that in contrast to the Eq. © the part of the CI © 
does not conserve the energy of localized quasiparticles 
and therefore provides the energy flow out of the vortex 
cores. 

Having in hand the expressions for the CI we can find 
the corrections to the distribution functions generated 
by the vortex motion. At first we find the anisotropic 
corrections taking into account only the localized part 
of the CI © St(f) = St loc (f), Our key idea is that to 
induce the non-equilibrium distribution in Gallilean in- 
variant system of localized quasiparticles the accelerated 
motion of the vortex is needed ii/0. The solution of the 
kinetic equations can be found in the usual way. Substi- 
tuting the ansatz / = / (0) (£) - (pu)^ + /W to the 
Eq.© we obtain the expression for the non-equilibrium 



correction: 



Tlo 



(j||P • u + 7_l[p x uj • z v J 



de 



(5) 



where z v is a unit vector determined by the direction of 
vorticity, ^ = 1/[1 + (w t ;oc ) 2 ] and j± = woTWYU- 

The force acting on the moving vortex from the en- 
semble of localized quasiparticles can be calculated sub- 



(3) stituting the solution ([5]) to the conventional exprcssio: 



,9. 



Fioc = -Afiiu - M±[z v x u] , 



where the masses are determined by 



Mi 



IK-L) 




Tlocl\\(± 



oe 



(6) 



(7) 



F.S. 



where TV is the particle concentration and the brackets 
denote the Fermi surface averaging. The first term in 
the Eq.© is determined by the inertial vortex mass My 
which in the limit t; oc = 00 coincides with the expres- 
sion obtained in collisionless regime^. This inertial term 
does not lead to the dissipation, that is to the trans- 
formation of the vortex kinetic energy to the heat. On 
the contrary the second term in the Eq.© can provide 
the dissipation. It is easy to see that for the circular 
vortex motion the period-averaged work of force Fi oc is 
non-zero due to the second term. In case of non-axially 
symmetric vortices the expressions for masses Afiin ) be- 
come more involved and take into account the interplay 
between bound fermions and the modes of rotational vor- 
tex core motion^. 

Let us consider the generic example of the single vortex 
line motion in the self-induced superfiuid velocity field de- 
scribed in terms of the propagation of Kelvin waves (KW) 
along the vortex line. KW play a key part in the relax- 
ation of superfiuid turbulence at low temperatures. The 
KW cascade characterized by the universal spectru m 18 i 20 
carries the turbulent energy from the main length scale 
of the vortex line separation to essentially lower length 
scales. 

In the local induction approximation (LIA) valid for 
the large Kelvin wavelength the motion of vortex line2i 
is determined by the linear tension energy, Magnus force 
Fm = — kN[z v x u] where k = 2i:h/m is a vorticity and 
the force © exerted on the vortex line by the localized 
excitations 

TO||U + m±[z v x U] + [z v x U] = v s TJ zz (8) 

where U(z,t) is the vortex displacement, v s = kA/Att 
and A = ln(£/Z c ) where l c is the cutoff parameter, toii,^ = 
A^iV (71W). The dispersion relation following from the 
Eq.© has the form 



ujk = ±k~i> s + (fc v s ) (im± — mil) 



(9) 



The above expression can be extrapolated to the smaller 
wavelengths 2 ^ by setting the cutoff parameter l c = 1/k. 
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Note that the corrections to the KW frequency of the 
same type as given by Eq.© are necessarily taken into 
account to break the integrability of LIA and launch the 
KW cascade (see e.g. Refl33). Thus the corrections 
given by Eq. ((9]) can in principle have significant influence 
on the theory of quantum turbulence. 

The expression ^ demonstrates that even in the ab- 
sence of the heat bath the non-trivial dynamics of local- 
ized excitation can provide the finite decrement which 
can be estimated as t^ 1 = mj_uij { ~ t,~* (w^/wo) 2 . The 
decrement depends on the temperature T; oc and can be 
finite even if the heat bath temperature is zero T^ei = 0. 

To calculate the temperature T\ oc we use the method 
analogous to the one employed in^ 2 - with the only dif- 
ference that the energy flow out of the vortex core is 
determined by the collisions described by the CI (@| and 
not by the collisions with phonons. Let us substitute the 
distribution function <j5j) into the Eq. ((T|) and take the 
average by quasiparticle momentum direction over the 
Fermi surface and time t. From the first term in the 
Eq.JT]) we obtain 

/ r i d f (1) \ d2 f i0) 

where a = (u x u y — u y ii x ) t /2 and jj_ = (j±)fs- We have 
neglected the term containing the time average (u 2 ) 4 as- 
suming the periodic vortex motion. Also wc multiply 
and integrate the Eq. (p} by energy. Making use of the 
energy conserving property of localized CI © we obtain 
the equation which determines the temperature T; oc 

/Ao /-Ao ^2 

eSti-di{f)de — aj± / e— T de (10) 
-A J-A de 2 

The obtained equation has the simple physical meaning 
expressing the balance between the heat produced by the 
accelerated vortex motion and the energy flow out of the 
vortex core carried by the emitted delocalized quasipar- 
ticles. 

The further analytical analysis can by carried out as- 
suming that Tioc <C Ao. By the order of magnitude we 
have ujq ~ t" 1 ~ Ao/(fcp£) therefore u!qTi oc <C 1 and 
7j_ ~ ujqTi oc . Taking into account the temperature de- 
pendence of the collision time t; oc we obtain the simpli- 
fied heat balance equation in the form 

a^/V c 3 = (A /T loc )e- A °/ T ^, (11) 

where V c = Ao/pp is the critical supcrfluid velocity. 

To estimate the value of the vortex core temperature 
Ti oc which can be realized in the conditions of quantum 
turbulence in 3 He-B we evaluate the factor a in (fTTj) as- 
suming that the vortex dynamics is regulated by the KW 
cascade. Introducing the complex field W — U x + iU y 

we obtain a = (iWW* + c.c\ /4. The time average 
can be expressed through the 'kelvon' occupational num- 
ber n k = L(\W k \ 2 ) where W = k~ 1 / 2 Y,k Wk(t)e ikz and 



L ~ l~ 2 is the vortex line length per unit volume, I is 
the vortex line separation. To obtain the estimation of 
a we can use the spectrum of kelvons nk = Atfc~ 3 (fc/)~ 7 
where 7 = 2/5 or 7 = 2/3 according to Refs. [l8| and 
I20I correspondingly. We have omitted the dimension- 
less prefactors which are not important in our deriva- 
tion. The above scaling law results in the estimation 
a ~ k k^(k c l)~^ where k c is the cutoff wave number of 
the KW cascade^ 3 .. The heat balance equation (ITT1) can 
be rewritten in the form 

(A /T ioc ) e -Ao/T !oo w (fcce)4 - 7( £ /0 ^ (12) 

Further we need to estimate the value of the cutoff 
wave number k c of the KW cascade. This can be ob- 
tained in the usual way by comparing the large-scale 
energy flow in a Kolmogorov cascade^ 3 . E t ~ k?1~ a to 
the value of the energy dissipated by the KW cascade 

Pkw = I 2 J k " T^sxdk where Ek — ^k^k- From the 
equation Pk w = E t we obtain the cutoff parameter 

kc - T 1 (e/0" (Ao/T ioc ) 2/J (13) 

where (3 = (4 — 7) _1 . Substituting the obtained value 
of k c to the (fT2]) and taking into account that for the 
typical experiments in 3 He-B £/l ~ 10 -4 we obtain the 
temperature Ti oc rj 0.1 Aq ~ 0.2T C . The value of the 
cutoff parameter is k c ~ 0.2^ _1 corresponding to the 
Kelvin wavelength A = 2-k /k c ~ 30^ which is much larger 
than the core size £ and therefore can be reached by the 
KW cascade. Thus we obtain that at temperatures lower 
than 0.2T C the dynamics of the vortices in 3 He-B should 
be determined by the local heating of the vortex cores 
rather than by the interaction with the heat bath. 

It is interesting also to compare the power of the en- 
ergy losses predicted by the theory with the experimen- 
tally observed values of the energy flux out of the tur- 
bulent region generated by the propagating vortex front 
in 3 He-B at the temperatures down to Q.2Tj^. We use 
the typical 3 He parameters £ = 10~ 6 cm, (&f£) ~ 10 2 , 
A ~ 10" 25 J,t„ ~ 10~ 6 s and assume T loc w 0.1 A to 
obtain that P kw ~ 10" 10 (fc c £) 4 "' ) 'VF/cm 3 . In the experi- 
ment the typical value of the turbulent region is ~ 1cm 
the heat flux is I E w 10~ 10 (fc c £) 4 - 7 W^/cm 2 which is close 
to the experimentally observed value if (fc c £) ~ 1. Note 
that this estimation is the minimal value of the heat flux 
which can be obtained in the limit T — > 0. 

In conclusion we have studied the dynamics of quan- 
tized vortices in ultra cold Fermi supcrfluids. The new 
mechanism is suggested of how the energy can be trans- 
ferred to the ensemble of localized quasiparticcs result- 
ing in the local temperature rise inside the vortex cores 
and dissipation due to the escape of hot quasiparticles 
into delocalized states. The force acting on vortex is de- 
rived and shown to provide the finite decrement of Kelvin 
waves which is determined by the vortex core tempera- 
ture. It would be interesting to investigate the dissipativc 
vortex dynamics under the action of this force in more 
detail in large-scale numerical simulations. 
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